Short-range exchange and correlation energy density functionals: beyond the local 

density approximation 



Julian Toulouse, Frangois Colonna, and Andreas SavirQ 
Laboratoire de Chimie Theorique, 
CNRS et Universite Pierre et Mane Curie, 
4 place Jussieu, 75252 Paris, France 
(Dated: February 2, 2008) 

We propose approximations which go beyond the local density approximation for the short-range 
exchange and correlation density functionals appearing in a multi-determinantal extension of the 
Kohn-Sham scheme. A first approximation consists in defining locally the range of the interaction 
in the correlation functional. Another approximation, more conventional, is based on a gradient 
expansion of the short-range exchange-correlation functional. Finally, we also test a short-range 
generalized-gradient approximation by extending the Perdew-Burke-Ernzerhof exchange-correlation 
functional to short-range interactions. 



I. INTRODUCTION 

The Kohn-Sham (KS) scheme pj of density functional 
theory (DFT) (see, e.g., 0,0) can be extended to han- 
dle multi-determinantal wave functions Q. This enables 
to describe near-degeneracy correlation effects by explic- 
itly taking into account several Slater determinants. The 
method is based on a decomposition of the true Coulomb 
electron-electron interaction as 



(1) 



where vg e (r) is a long-range interaction and vg e (r) is the 
complement short-range interaction. This separation is 
controlled by the parameter /i. In previous studies 5, 6, 
0, the error function has been used to define the long- 
range interaction 



erf (fir) 



(2) 



referred to as the erf interaction. More recently [? ] , we 
have proposed a sharper long-range/short-range separa- 
tion with the erfgau interaction 



v 



ee, erfgau 



(r) = 



erf(c^r) 2c/i 



(3) 



where c = (l + 6^3) « 3.375 is a constant chosen to 
facilitate the comparison with the erf interaction. In- 
deed, with this choice, the parameter \i has roughly the 
same physical meaning for the erf and erfgau interac- 
tions: represents the distance beyond which the in- 
teraction reduces to the usual Coulomb long-range tail. 
(Note that in Ref. ? the constant c was absorbed 
into the parameter fi). Both modified interactions de- 
fine a generalized adiabatic connection Q between the 
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non-interacting KS system corresponding to fi = and 
the physical interacting system corresponding to \± — ► oo. 

The ground-state energy of an electronic system can 
then be (exactly in principle) calculated by the sum of a 
wave function part associated to the long-range electron- 
electron interaction V£, = J2i<j v eX r ij) an( ^ a density 
functional part describing the remaining interactions 

E = (*"|T + V£ + V W \W) + U»[n] + E£ c [n], (4) 

where T is the kinetic energy operator, V ne = iw( r i) 
is the nuclei-electron interaction, is the short-range 
Hartree energy and E£ c is the short-range exchange- 
correlation functional defined as the difference between 
the standard KS exchange-correlation energy E xc and the 
long-range exchange-correlation energy E£ c associated to 
the interaction v£. 



E^ . 



(5) 



In Eq. QJ, ^ is the multi-determinantal ground-state 
wave function of a fictitious system with long-range 
electron-electron interaction VJt 



f + VJi + Vn =£^1^) 



(6) 



where = ^ M ( r i) is the external local potential en- 
suring that this fictitious system has the same ground 
state density n than the physical system. The potential 
is simply related to the functional derivative of the 
short-range functionals through 



X*) 



(5l> 5E» 



6n(r) Sn(r) 



(7) 



Previous applications of the method [y, |Jj show that, 
for a reasonable long-range/short-range separation (/z not 
too small) and if the few most important (nearly degener- 
ate) configurations are included in the wave function 
good results are obtained for total atomic and molecu- 
lar energies, including near-degenerate systems. This is 
remarkable since, in these previous studies, only the sim- 
ple local density approximation (LDA) was used for the 
short-range functional E£ c [n]. 
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The purpose of the present work is to further improve 
the accuracy of the method by proposing approxima- 
tions for Eg c [ri\ which go beyond the LDA. We will not 
present in this paper total energies E (or equivalently 
total exchange-correlation energies E xc ), but will focus 
only on short-range exchange-correlation energies E£ c . 
Indeed, for a chosen /i, it is the approximation for Eg c [n] 
that limits in practice the accuracy of the method, not 
the other long-range contributions to the energy that can 
be calculated by a traditional ab initio wave function 
method. We will test our proposed approximations for 
Eg c by comparison to accurate data obtained for small 
atomic systems. 

The paper is organized as follows. In Sec. |nj we dis- 
cuss the limitations of the LDA for Eg c [n]. In Sec. II I II 
we introduce a local interaction parameter fj, to improve 
the LDA correlation functional. We then discuss the ex- 
tension of traditional gradient corrections to the short- 
range functional E!£ c [n]. Indeed, the short-range versions 
of the gradient expansion approximation (GEA) is pre- 
sented in Sec. lIVI while Sec. Ivl contains the extension of a 
generalized-gradient approximation (GGA), namely the 
PBE functional :Q], to a modified interaction. Finally, 
Sec. IVII contains our concluding remarks. 

Atomic units will be used throughout this work. 



II. LIMITATIONS OF THE LDA 

In the LDA, the short-range exchange-correlation en- 
ergy per particle £^ mif W of a uniform electron gas with 
modified interaction 0, is transferred locally to the 
inhomogeneous system of interest 



Eg 



n(v)^r if (n(r))dr. 



(8) 



To underline the dependence on the interaction parame- 
ter, we will refer to this approximation as /i-LDA. 

In Figs. Hand [3 we have plotted the short-range LDA 
exchange and correlation energies of the Be atom with 
respect to /i for the erf and erfgau interactions. An ac- 
curate density obtained from a multi-reference configura- 
tion interaction calculation with single and double excita- 
tions (MRCISD) [H[ll 

has been used. For comparison, 
accurate calculations of the exchange and correlations en- 
ergies along the adiabatic connections are also reported. 
In these accurate calculations, we start from an accurate 
reference density and, for each /i, numerically optimize 
the external potential u^(r) appearing in Eq. © so as 
to recover the reference density. The ground-state wave 
function vf^ is then computed according to Eq. © by 
MRCISD and the various energy components like the 
short-range exchange and correlation ener gies are de- 
duced. (For further details, see 

Refs.mElllil). 

For both the erf and erfgau interactions, the /i-LDA is 

very accurate for large /j, but fails near the KS end (// = 0) 
of the adiabatic connection. In particular, the exchange 
energy is underestimated and the correlation energy is 
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FIG. 1: Accurate (full curves) and /i-LDA (dashed curves) 
short-range exchange energy along the erf (thick curves) and 
erfgau (thin curves) adiabatic connections for the Be atom. 
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FIG. 2: Accurate (full curves) and /i-LDA (dashed curves) 
short-range correlation energy along the erf (thick curves) 
and erfgau (thin curves) adiabatic connections for the Be 
atom. 



overestimated. Thanks to a sharper separation of long- 
range and short-range electron-electron interactions, the 
erfgau interaction provides a slight improvement over the 
erf interaction. In fact, Fig. shows that the erfgau \x- 
LDA exchange energy curve reaches the exact one for a 
larger energy (ss —0.75 Hartree) and a smaller /i (w 2.5) 
than in the erf case (ss —0.5 Hartree at /i ss 3), meaning 
that the erfgau /i-LDA is able to correctly describe a 
larger part of the exchange energy. For the correlation 
energy, the two interactions lead to similar results. We 
note in passing that the unimportant little bump on the 
curve of the accurate correlation energy for small /i with 
the erfgau interaction in Fig. [21 (and in other following 
figures) is a manifestation of the non-monotonicity of this 
interaction with respect to /i (being in turn responsible 
for its attractive character for very small /i [13 )• I n the 
remaining of the paper, we will present results only for 
the erfgau interaction. 

The performance of the /i-LDA can be further ana- 
lyzed by looking, for instance, at the correlation energy 
density s£(r). There is not unique definition of this quan- 
tity; we used the definition of Ref. or based on the 
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integration over the adiabatic connection 

£t(ri) = \ f <H J dran^n.ra)^^, (9) 

where n|(ri,r2) is the correlation hole for the interac- 
tion parameter £. In Fig. |3| this correlation energy den- 
sity accurately computed for the Be atom is compared 
with the /i-LDA with the erfgau interaction for a series 
of ju's. For fi = (KS system), the correlation energy 
density of Eq. © is largely overestimated by the LDA 
over the whole range of r. When /i is increased (Fig. [3| 
with fj, — 0.21 and /j, = 1.20), the /(-LDA starts to bet- 
ter reproduces the accurate energy density in the valence 
region (r > 1) but still overestimates it in the core re- 
gion (r < 1). For the /i-LDA to also well reproduce the 
core region, larger values of \i are required (Fig. [3] with 
M = 3.00). 

Having in mind that the choice of a large value of /i 
means an important computational effort for the part of 
the energy treated by wave function methods, it is im- 
portant to improve the short-range exchange-correlation 
functional toward the small /i region of the adiabatic con- 
nection. 



III. LDA CORRELATION WITH LOCAL 
INTERACTION PARAMETER 

In inhomogeneous systems the electrons are correlated 
on a shorter distance than in the uniform electron gas. 
Thus, within the KS scheme, the LDA transfers spurious 
long-range correlation contributions from the electron gas 
to the finite system of interest. This point have long been 
understood and have guided important developments in 
density functional approximations. For instance, Lan- 
greth and Mehl 0,0] proposed a gradient correction to 
the LDA based on a cut-off in momentum-space for small 
k which basically corresponds to removing the long-range 
interaction effects. 

In the method proposed here, for small /i, i.e. when a 
large range of the electron-electron interaction is retained 
in the functional part, the /i-LDA also transfers spurious 
long-range interaction effects from the uniform electron 
gas to finite systems and leads to an overestimation of the 
correlation energy. A simple way to eliminate this incor- 
rect long-range interaction effects, is to use in the /i-LDA 
correlation functional an effective interaction parameter, 
larger than the one used in the wave function part of 
the calculation. Of course, in the uniform electron gas, 
the range of the interaction is relative to the density (the 
relevant scaled variable is r/r s with r s — (47m/3)~ 1//3 ) 
and consequently one has to actually consider an effec- 
tive local interaction parameter, /iz(r), defined from the 
inhomogeneous density profile n(r). However, the in- 
teraction parameter used in the functional must not be 
smaller that the global \x used in the wave function part 
of the calculation to avoid a double counting of correla- 
tion effects. We therefore take in the correlation /i-LDA 



functional an effective local interaction parameter given 
by 

/ieff(r) = max(/i;(r ),/*). (10) 

Thus, this approximation corresponds to simply drop off 
some long-range correlation from the LDA in regions of 
space where /i;(r) > /i. 

Several choices are possible for the local interaction 
parameter /i;(r). The first natural length scale in real 
space for the electron-electron interaction is provided by 
the Wigner-Seitz radius r s , representing the radius of the 
sphere centered in r containing one electron. Equiva- 
lently, the same length scale can be expressed in momen- 
tum space by the Fermi wave vector k F = 1/ (ar s ) where 
a = (97r/4) -1 / 3 . Following a previous work pfl, we thus 
take as a local interaction parameter 

W,a(r) = —rx = ak F (r). (11) 
r s {r) 

However, studies of the uniform electron gas indicate that 
the relevant length scale for correlation effects in momen- 
tum space is given by the Thomas-Fermi wave vector 
k s = y 'Akpj-K describing the screening of the Coulomb 
interaction. Similarly to Eq. 1)11(1 . we can thus take as 
local interaction parameter 

W>b (r)=aA.(r) = H^* (12) 

More elaborated dependences of [ii on r s are possible. 
For example, Ferrell |l9j | has determined a momentum 
cut-off in the uniform electron gas corresponding to the 
shortest wavelength at which plasma oscillations can be 
sustained by the electron gas. It turns out however that 
results obtained with a local interaction parameter set 
to Ferrel's momentum cut-off are very similar to those 
obtained with the local interaction parameter of Eq. 1(12(1 
and will not be discussed here. 

A much more interesting possibility is to choose a lo- 
cal interaction parameter depending on the gradient of 
the density. In fact, as summarized by Svendsen and von 
Barth HJ from the work of Langreth and Mehl [H 
and Perdew et al. 2\\, "distances over which electrons 
are correlated in realistic inhomogeneous systems are de- 
termined more by actual variations in the density profile 
than by the screening of an electron gas at the local value 
of the density". We use a simple geometrical argument 
similar of that of Langreth and Mehl to determine a 
local characteristic distance d(r) over which the density 
varies in an inhomogeneous system. Consider a density 
distribution with constant density gradient consisting of 
an isosceles triangle as a prototype of an atom. The ex- 
tension of this density distribution is d = 4n/|Vn| where 
n is the density at the middle of one side of the trian- 
gle, i.e. a typical value of the density. We thus define a 
(semi)local interaction parameter by 

|Vn(r)| 

w ' c(r) = "toST (13) 
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FIG. 3: Short-range correlation energy density s%(r) with respect to r for the Be atom with the erfgau interaction for /i = 0, 
0.21, 1.20, 3.00. Accurate data computed from integration over the adiabatic connection (Eq.|5| full curve) are compared with 
the /i-LDA (dashed curve). 



The constant 4 in Eq. (|13f) is somehow arbitrary; experi- 
ence shows that this is reasonable. 

Fig.0]shows the correlation energy is greatly improved 
by the local [i approach. However, for the choices /i; )a = 
akp or /x/.b = ctk Sl a more detailed analysis through 
the correlation energy density in Fig. [3] indicates the [i- 
LDA functional with a local fi is only an average of the 
"exact" correlation energy density. Of course, the /i-LDA 
functional with a local fj, improves on increasing of fj, (not 
shown) but it is clear that both choices are inadequate 
for describing the shell structure. On the contrary, one 
sees on the same Fig. that the choice /// iC = |Vn|/(4n) 
enables to recover well the shell structure in Be. We 
therefore consider this last choice as the more appropriate 
local parameter fii for improving the LDA correlation 
functional. 

Finally, we note that this local fj, approach cannot be 
directly applied to the exchange energy. Indeed, the LDA 
underestimates the exchange energy and thus choosing a 
large effective interaction parameter in the /i-LDA would 
only deteriorate the exchange energy. 



IV. GRADIENT EXPANSIONS 

We now consider more conventional approaches to cor- 
rect the local density approximation at small interaction 
parameter, namely gradient corrections. In the gradi- 
ent expansion approximation (GEA), the short-range ex- 
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FIG. 4: Short-range correlation energy along the erfgau adia- 
batic connection for the Be atom. Accurate data (full curve) 
are compared to the /i-LDA functional with a global fi (thick 
long-dashed curve) and with (semi) local interaction param- 
eters: /Ltj >a = otkF (long-dashed curve), /i;.b = ak a (short- 
dashed curve) and fj,i tC — Vn|/(4n) (dotted curve). 



change energy density is expanded in term of the reduced 
gradient s = |Vn|/(2fej?n) and the reduced Laplacian 
q = V 2 n/{{2k F ) 2 n) 

e£' GEA - Lap (n) = e^ unif (n) [l + a(/2)s 2 + d(fi)q] , (14) 

where the reduced interaction parameter /i = fi/(2kp) 
has been introduced. The gradient and Laplacian coef- 
ficients a(fl) and d(fl) can be calculated numerically for 
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FIG. 5: Short-range correlation energy density e^{r) with 
respect to r for the Be atom with the erfgau interaction 
for /i = 0.2f . Accurate data (Eq. |U] full curve) computed 
from integration over the adiabatic connection are compared 
with the /i-LDA functional with local interaction parameters: 
Mi, a = akF (long-dashed curve), /i;.b = atk s (short-dashed 
curve) and /ii, c = Vn|/(4n) (dotted curve). 

all /i's except for /i = where there are not defined (see 
Appendix . Alternatively, one can perform an inte- 
gration by parts and obtain the more conventional GEA 
exchange density 

^ GEA (n)=er nif (n)[l + &(A> 2 ], (15) 

where the exchange second-order gradient coefficient b(fi) 
has been calculated analytically for the erf and erfgau in- 
teractions [see Appendix^] Eqs. l|A18f) and i|A19fl ]. The 
approximations of Eqs. Q14f) and (|15f) will be referred to 
as /i-GEA-Lap and /i-GEA, respectively. 

Fig. compares the /i-GEA-Lap and /z-GEA short- 
range exchange energy densities for the Be atom with 
the erfgau interaction at /z = 0.21. The exchange energy 
density in the /i-LDA and calculated from the accurate 
exchange hole n x (ri,r2), 

et(ri) = \J rfr 2 n x (r 1; r 2 X(r 12 ), (16) 

are also reported. To avoid the Laplacian divergence 
at the nucleus in /i-GEA-Lap, radial energy densities 
A r Kr 2, n{r)ef.(r) are actually represented. One sees that 
the integration by parts in /z-GEA does not change qual- 
itatively the shape of the curve and it is meaningful to 
compare both the /i-GEA-Lap and /i-GEA exchange en- 
ergy densities to the accurate and /i-LDA ones. Besides, 
for larger values of /z, /i-GEA-Lap and /i-GEA become 
nearly identical (not shown). Of course, /i-GEA-Lap and 
/i-GEA give the same short-range exchange energy 
and we will used in practice /i-GEA because of its great- 
est simplicity. 

Similar to Eq. (|15|l . the short-range correlation energy 
density is expanded in term of the reduced gradient t = 
\Vn\/(2k s n) 

s^ GEA (n) = eZ> unil (n) + /3(r 5 ,/z)i 2 . (17) 



FIG. 6: Radial short-range exchange energy density 
47rr 2 n(r)e^(r) with respect to r for the Be atom with the 
erfgau interaction for /i = 0.21. Accurate data fEq. 1161 full 
curve) are compared with the /i-LDA (long-dashed curve), 
the /i-GEA (short-dashed curve) and the /i-GEA-Lap (dotted 
curve) . 

The correlation second-order gradient coefficient f3(r s ,fi) 
can be estimated by imposing that the correlation gra- 
dient expansion cancels the exchange gradient expansion 

[3(r s , M ) = -e^{n)b{fi){ks/k F f. (18) 

This is based on the observation that, for /i = 0, the 
exchange-correlation contribution to the linear response 
of the uniform electron gas to a slowly oscillating pertur- 
bation (small momenta) almost vanishes • This short- 
range GEA functional approximation will be referred to 
as /i-GEA. Our interest is to evaluate the improvement 
brought by /i-GEA over /i-LDA for the exchange and cor- 
relation energies separately. Notice that, in our construc- 
tion, if we consider exchange and correlation together /i- 
GEA is identical to /i-LDA. 

The /i-GEA exchange energy is plotted in Fig El F° r 
large /i (/i > 3), the range of the interaction is too small 
to feel the slowly oscillating gradient correction and the 
/i-GEA functional reduces to the /i-LDA functional. The 
/i-GEA curve reaches the exact one at a smaller value of 
fi (/i w 2) than the /i-LDA curve (/z ~ 3). For small /i, 
the /i-GEA reduces the /i-LDA error by about a factor 
two. For the correlation energy (Fig.[8J), the /i-GEA also 
extends the domain of validity of the /i-LDA when /i is 
decreased down to /z ~ 4. For smaller /i when long-range 
correlation effects are introduced in the functional, the 
gradient expansion breaks down. 

Whereas the GEA for the KS scheme (/i = 0) often 
disappointingly deteriorates the LDA in real inhomoge- 
neous systems [23, 01 > ^ can be consider as a systematic 
way to extend the domain of validity of the /i-LDA for 
the short-range part only of the exchange and correlation 
functionals. These results are in agreement with the wave 
vector analysis of the GEA by Langreth and Perdew j2f| 
who found that the GEA works well for large momenta 
(corresponding to short-range density fluctuations). 
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FIG. 7: Short-range exchange energy along the erfgau adia- 
batic connection for the Be atom. Accurate data (full curve) 
are compared to the /i-LDA functional (long-dashed curve), 
the /i-GEA functional (dotted curve) and the /i-PBE func- 
tional (short-dashed curve) . The curves corresponding to the 
/i-GEA and /i-PBE functionals are nearly superimposed. 
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FIG. 8: Short-range correlation energy along the erfgau adia- 
batic connection for the Be atom. Accurate data (full curve) 
are compared to the /i-LDA functional (long-dashed curve), 
the /x-GEA functional (dotted curve) and the /i-PBE func- 
tional (short-dashed curve) . 



V. SHORT-RANGE PBE FUNCTIONAL 

Starting from the gradient expansions of the previous 
section, it is possible to construct generalized-gradient 
approximations (GGA) for the short-range exchange and 
correlation energy functionals which cure the main short- 
comings of the /i-GEA. We have actually constructed a 
modified version of the PBE functional 9] for the mod- 
ified interactions. All of the details of the construction 
of this functional that we will call /i-PBE are given in 
Appendix [BJ We simply mention here that we have use 
the same ansatz than PBE and impose similar theoreti- 
cal constraints. However, a difference with the PBE ap- 
proach is that we impose the second-order gradient coef- 
ficient for exchange and approximate that of correlation 
from it while PBE did the opposite. Consequently, our 
/i-PBE functional does not reduce to the original PBE 
for /i = 0. 
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FIG. 9: Short-range correlation energy density e£ (r) with re- 
spect to r for the Be atom with the erfgau interaction for 
/i = 0.21. Accurate data (full curve) computed from inte- 
gration over the adiabatic connection are compared with the 
/i-PBE functional (dashed curve). 



The exchange and correlation energies of the Be atom 
obtained with this /i-PBE functional along the erfgau 
adiabatic connection are compared to the /i-LDA and /i- 
GEA functionals in Figs. and [S] For a simple system 
like the Be atom, the /i-GEA and /i-PBE functionals are 
nearly identical for the exchange energy. For the corre- 
lation part, one sees that /i-PBE greatly improved the 
/i-LDA and /i-GEA results. The /i-PBE correlation en- 
ergy is indeed very accurate along the entire adiabatic 
connection, except near the KS end (/i = 0) where a 
significant error remains. This inaccuracy of (semi) local 
functionals at /i = is in favor of the strategy consisting 
in escaping the KS scheme by using a short-range func- 
tional at finite /i and treating the remaining part of the 
energy by other, more appropriate methods. In Fig. [5J 
it is shown that the correlation energy density of /i-PBE 
functional qualitatively describes the shell structure of 
the Be atom for small /i. 

Finally, the difficult case of the Ne 6+ atom, presenting 
strong near-degeneracy correlation effects, is reported in 
Fig- El The correlation energy given by the /i-PBE func- 
tional actually over-corrects that of the /i— LDA. Again, 
this result suggests that for small \x semilocal approxi- 
mations like the /i-PBE functional are inappropriate, es- 
pecially when near-degeneracy correlation effects play an 
important role. 



VI. CONCLUDING REMARKS 

In this paper, we have considered exchange and cor- 
relation energy functionals associated to short-range 
electron-electron interactions. We have proposed new 
approximations which partially correct the local density 
approximation. Indeed, we showed that the short-range 
LDA correlation energy can be significantly corrected by 
simply adjusting in the functional the range of the in- 
teraction locally from the density and its gradient. We 
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FIG. 10: Short-range correlation energy along the erfgau adi- 
abatic connection for the Ne 6+ atom. Accurate data (full 
curve) are compared to the fi-LDA functional (long-dashed 
curve) and the /i-PBE functional (short-dashed curve). 

have also constructed generalizations of the GEA and 
of a GGA (namely, the PBE functional) to the case of 
short-range interactions. These new short-range func- 
tional treat well a larger range of electron-electron inter- 
action than the short-range LDA. However, better (possi- 
bly non-local) approximations are still needed to further 
extend the range of interaction well treated by the func- 
tional. 



APPENDIX A: GRADIENT EXPANSION OF 
THE SHORT-RANGE EXCHANGE ENERGY 
FUNCTIONAL 

The short-range exchange energy functional writes 

E£ = ~ / 1 n(r)n x (r,r + R)v£ e (R)drdR, (Al) 



The exchange energy density is 

£ -GEA-Lap Ks ^ ~ } = 1 J n GEA (l . r + R) -p {R)dR 
If 00 Z 

(A12) 



where n x (r, r + R) is the exchange hole. In the gradient- 
expansion approximation (GEA), the exchange hole is 
expanded to second order in term of the gradients of the 
density pll^ 

«f A (r, r + R) = ~n(v)y(r, r + R), (A2) 

where 

y = J + Lk F 3 R-Vk 2 F + Mk F 6 (K-\7k 2 F ) 2 
+Nk F 6 (Wk 2 F ) 2 + zLk F 4 (R • W) 2 k 2 F /6 
-z 2 Jk F 4 V 2 k 2 F /A8, (A3) 

with k F (r) = (37r 2 n(r)) 1 / 3 , R = R/R, z = 2k F R and 
J{z) = 72[4 + z 2 - (4- z 2 ) cos z-4z sin z]/z 6 , (A4) 



L(z) = 9(2-2cosz-zsinz)/z 3 , (A5) 



M(z) = 9(-z cos z + sin z)/(16z), (A6) 



N(z) = 3[8-(8-4z 2 )cosz-(8z-z 3 )sinz]/(16z 4 ). (A7) 

The spherical average of the GEA exchange hole in term 
of the reduced gradient s — \Vn\ / (2k F n) and the reduced 
Laplacian q — V 2 n/ ((2k F ) 2 n) writes 



(A8) 



(A9) 



(A10) 



(All) 



where fi — fi/(2k F ). For convenience, e^? EA " Lap (n, s, q, fx) 
can be expressed by 

£ -GEA-Lap (nj fi] 9j -J = £ -unif (n; ~) ^ + ^2 + _ 

(A13) 



n^ EA (k F ,s,q,z) = ±-J nf A (r,r + R)dSl R = n^ a (k F ,z) + n^ ad (k F , z)s 2 +n l ^(k F ,z)q, 



with the exchange hole of the uniform electron gas 



if 24fc 3 [zcos(z/2)-2sin(z/2)] 2 

U x [K Fl Z)— -jr- - . 



the gradient coefficient 



d _ fc|[-4(18 + z 2 ) + (72 - 32z 2 + 3z 4 ) cosz + 2z(36 - 5z 2 ) sinz] 

and the Laplacian coefficient 

fcf,[12(12 + z 2 ) + 12(-12 + 5z 2 ) cosz + 12z(-12 + z 2 ) sin z\ 



4 ap (fc F ,z) 



54tt 2 z 4 
I 
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where e" nlf (n, fi) has been given in RefeJllIS (Note 
that, for the erfgau interaction, in Ref. [10| the constant 
c of Eq. (PJ was absorbed into the parameter /i). The 
integral in Eq. (|A12|) diverges for /i = due to the 
long-range character of the Coulomb interaction and thus 
£^f EA ~ Lap (n, s,q,jl = 0) is not defined. (For a recent dis- 
cussion on the non-analyticity of the inhomogcncous ex- 
change energy density with respect to the electron den- 
sity, see Refs. |2^, |29|) . Nevertheless, for finite /i, the 
integral does exist and the gradient and Laplacian coef- 



ficients, a(fl) and d(fi), can be computed numerically. 

To eliminate the divergence at fi = 0, one can per- 
form an integration by parts of (k F , s, q, z) over r 
in Eq. (|A1|) and define a new GEA exchange hole in term 
of the reduced gradient s only 

n^ EA (k F , s, z) = < nif (fc F , z) + nf ad (k F , z)s 2 , (A14) 

with the associated gradient coefficient 



~grad 



(k F ,z) 



k 3 F [-72 + (72 - 36z 2 + z 4 ) cos z - 2z(-36 + 5z 2 ) sin z] 
54tt 2 z 4 : 



(A15) 



and the corresponding exchange energy density 



e^ A (n,s,fi) 
1 



2(2k F y 



(k F , s, z)v% e {-—)A-K 2 z 2 dz. 
2k F 



(A16) 

Similarly to Eq. (|A13|I . e^ EA (n, s, fi) is expressed as 

e^ A (n, a, fi) = ^(n, ft [l + b{ft)s 2 ] . (A17) 

The integral in Eq. (|A16I) is now defined for all /i's and 
can be done analytically for both the erf and erfgau mod- 
ified interactions. For the gradient coefficient b(p.) with 
the erf interaction, we obtain 



&erf(£) 



-Ci + c 2 e 



l/(4/I 2 ) 



c 3 + hAaeWP) ' 
and, with the erfgau interaction, 

-3V3c 5 + 2^ 2 c 1 e 1 /(2- 2 ) + C6 e 3 /( 4 ^ 2 ) 



(A18) 



^erfgau (A) 



c 7 - 2i/ 2 c 3 e 1 /(2^ 2 ) + 12i/ 2 (-9c 4 + c 8 )e 3 /( 4 ^ 2 ) ' 

(A19) 

with ci = 1 + 22^ 2 + 144^ 4 , c 2 = 2^ 2 (-7 + 72v 2 ), 
c 3 = -864i/ 4 (-l + 2^ 2 ), c 4 = ^ 2 [-3 - 24i/ 2 + 32^ 2 + 
8iVierf(l/(2i/))], c 5 = 3 + 18^ 2 + 48^ 4 + 64z/ 6 , c 6 = 
4i/ 4 (7 - 72// 2 + 48Vli/ 2 ), c 7 = -192V3^ 6 (-3 + 8z/ 2 ), 
! [-18V3^+16\/3^ 3 +9 v ^ : erf(V3/(2iy))], v = fi/c 



cs 



and c = (l + 6\/3) ~ 3.375. For // = for both inter- 
actions, b(fj) reduces to Sham's coefficient [3(| El of the 
standard exchange GEA 
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81' 



6(/i = 0) = — . (A20) 

The approximations e^ EA " Lap (n, s, g, /2) and 
(n, s, fi) will be referred to as /i-GEA-Lap and 



=GEA 



/i-GEA, respectively. They are compared in Sec. hV\ 
for the case of the Be atom with the erfgau interaction 




7.5 10 12.5 15 17.5 



FIG. 11: z-integrand of the short-range /i-GEA exchange en- 
ergy density e^ EA (n, s,/i) fEo, IA16t with the erf interaction 
and with s — and /t — (full curve), s = 1 and /i = (long- 
dashed curve) and s = 1 and \i = 0.2 (short-dashed curve), 
all for kF = 1. 



(Fig. 0) . Because of its greatest simplicity, 



.-GEA 



[n,s,p) 



will be used in practice. 

To see the effect the short-range interaction v£ e (R), 
the z-integrand of the /i-GEA exchange energy density 
(Eq. IA16J) has been plotted with respect to z in Fig. El 
in the case of the erf interaction. One sees that the 
spurious strong oscillations of the GEA exchange hole 
at large interelectronic distances are efficiently cut off at 
finite \i by the short-range interaction. In other words, 
only the short-range part of the GEA exchange hole is 
used in the short-range exchange functional. 

In order to appreciate the contribution of the gradient 
term of the /i-GEA with respect to the interaction pa- 
rameter, we have represented e™ (n,/2)6(/i)/e" nlf (n) in 
Fig. 1121 for the erf and erfgau interactions. Not surpris- 
ingly, the gradient correction is decreased for a short- 
range interaction. Notice however that for small /2, the 
gradient correction for the erfgau is increased compared 
to the Coulomb case fx = 0. 
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0.4 0.6 

fj/2k F 

FIG. 12: Ratio of the gradient contribution of /i-GEA short- 
range exchange energy density to the LDA exchange energy 
density with Coulomb interaction e" nlf (n, ji)b(jl) / e" 111 (n) with 
respect to /z/(2fcfr) for the erf (full curve) and erfgau (dashed 
curve) interactions. The horizontal line corresponds to the 
Coulomb case (/z = 0) where fe(/i = 0) = 7/81. 



APPENDIX B: SHORT-RANGE PBE 
FUNCTIONAL 



In the context of Kohn-Sham DFT, Perdew, Burke and 
Ernzerhof Q have proposed a gradient-corrected func- 
tional where the parameters are determined so as to sat- 
isfy theoretical conditions. We propose here to gener- 
alize the PBE functional along the erf or erfgau adia- 
batic connection. We note that for a different purpose 
Heyd, Scuseria and Ernzerhof [33| has already derived a 
short-range version of the exchange PBE functional cor- 
responding to the erf interaction by using a model of 
the PBE exchange hole |33|. (The obtained functional is 
called wPBE where u> corresponds to our parameter /z). 
We shall not follow this procedure but come back instead 
to the original formulation of the PBE functional. 

We thus take for the short-range exchange energy func- 
tional the PBE form 



E£ = J drn(r) £ f f (n(r) ) /i(r))F I ( g (r) ) /i(r)) ) (Bl) 

where we have introduced a dependence on the reduced 
interaction parameter p, — Li/{2kp). In Eq. IjBll) . s — 
|Vn|/(2fcpn) is a reduced density gradient, e" mf (n,/i) is 
the exchange energy per particle of a uniform electron 
gas with short-range interaction (see Refs. 0, Hoi) , and 
F£{s,p) is the enhancement factor 



F x (s,p) = 1 + k(P) 



k(/z) 



l + b(fl)s 2 /K(fj,) 



(B2) 



Note that this form had also be proposed by Becke |34j . 
The constant b(p.) is fixed by imposing the correct gradi- 
ent expansion of F x (s, fx) as s — > 



o . 1 



0.2 0.3 

^/2k F 



. 4 



0.5 



FIG. 13: Coefficients n of the /i-PBE exchange functional 
with respect to fj,/(2kp) for the erf (full curve) and erfgau 
(dashed curve) interactions. 



while k(/2) is connected to the rapidly varying limit (s 

oo) 



F x (s,fi) -> 1 + k(Ji). 



(B4) 



For b(fi) we take the coefficient arising from the gradi- 
ent expansion of the short-range exchange energy with 
erf or erfgau interaction (see Appendix [S] Eqs. IA18I 
and IA19|) . k(/2) is determined by imposing the Lieb- 
Oxford bound [35( which still holds for the short-range 
exchange functional E£ 



Eit>E x > -C 



t(r) 4 / 3 dr, 



(B5) 



since E£ = E x — E^. and the long-range exchange energy 
E^ is always negative. The constant C, for which Lieb 
and Oxford originally found 1.6787, has recently be im- 
proved by Chan and Handy to a value C — 1.6358 [3(|. 
A sufficient (but not necessary) condition for this bound 
to be satisfied is 



F s ( S ,M)<-CnVVerW). 



(B6) 



We thus take nQi) = -Cn 1/3 /e™ if (n. A) ~ !> th e largest 
value insuring condition (|B6|I . Fig. ll3l shows k(l1) for the 
erf and erfgau interaction. One sees that k(/2) increases 
with jl, but the Lieb-Oxford bound is actually of no effect 
for large fi since the GEA gradient correction vanishes 
anyway and the enhancement factor reduces to 1. 
The short-range correlation energy is given by 



EH 



drn(r)[ £ -r t (r s (r),M)+^(r s (r),f(r), M )], (B7) 



with the reduced density gradient t = \Vn\/(2k s n). We 
imposed the same conditions on the gradient correction 
H(r s ,t,Li) as PBE. 

(a) The second-order gradient expansion in the slowly 
varying limit (t — > 0) is 



F x (s,fi)~l + b{fl)s 2 + 



(B3) 



H{r„t,n) ^p(r s ,^)t 2 + 



(B8) 
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where the coefficient P(r s ,fi) is estimated from the ex- 
change gradient expansion according to Eq. lfT5)l. 

(b) The correlation is set to zero in the rapidly varying 
limit (t — > oo), thus 



H(r 3 ,t,fi) 



-unif 



(r s ,/i). 



(B9) 



(c) Under uniform scaling to the high-density limit, 



n(r) — » A J n(Ar), r s 



A 1/2 i and 



i.e 

A — > oo, the correlation energy scales to a constant. 
This well-known property for the Coulomb interaction 
case [33 is also valid for the short-range correlation 
functional. Thus, H(r s ,t,fi) must cancel the logarith- 
mic divergence of the local short-range correlation en- 
ergy e" nif (A _1 r s ,/i) « 7 ln(A- 1 ) + • ■ ■ with 7 w 0.031091, 
appearing as in the Coulomb case. Thus, as A — > 00, 



ff(A- 1 r s ,A 1 / 2 t,M)« 7 ln(A) + -- - . 
We then take the same ansatz than PBE 

P(r s ,[i)t 2 ( 1 + At 2 



(BIO) 



H(r s ,t,fx) =7 In 



1 



7 



1 + At 2 + A 2 t 4 



(BH) 



with 



1 



7 



df (r s ,M)/-y - 1 ' 



(B12) 



Note that in the original PBE functional, the second- 
order gradient coefficient for correlation was imposed in 
the high-density limit, (3(r s ->■ 0, /j, = 0) = j3 = 0.066725 
and the second-order gradient coefficient for exchange 
was chosen so that to exactly cancel the correlation gra- 
dient expansion: b(p, = 0) = /?(tt 2 /3) = 0.21951. On the 
contrary, we have used the exact exchange gradient ex- 
pansion for finite /i and estimate the correlation gradient 
expansion from it. Therefore, our //-PBE functional does 
not reduce to the original PBE functional for fi = 0. 
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